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INTRODUCTION 



The Landau theory of Fermi liquids has a long his- 
tory and remarkable results in describing the properties 
of electron liquid in ordinary metals and Fermi liquids of 
■^He type. The theory is based on the Landau paradigm 
that elementary excitations determine the physics at low 
temperatures. These excitations behave as quasiparti- 
cles, have a certain effective mass M*, which is indepen- 
dent of temperature T, number density x, and magnetic 
field strength B and is a parameter of the theory 
The discovery of strongly correlated Fermi systems rep- 
resented by heavy-fermion (HF) metals and 2D •^He ex- 
hibiting the non-Fermi liquid (NFL) behavior has opened 
tremendous challenges in the modern condensed matter 
physics [2|, 3., A, II, lS] ■ Facts collected on HF metals and 
2D ^He demonstrate that the effective mass strongly de- 
pends on T, X, B etc, while M* itself can reach very 
high values or even diverge 0, Such a behavior is 
so unusual that the traditional Landau quasiparticles 
paradigm does not apply to it. 

There is a common wisdom that quantum criticality, 
describing the collective fluctuations of matter under- 
going a second-order phase transition at zero temper- 
ature, suppresses quasiparticles and thus generates the 
NFL behavior, depending on the initialground state, ei- 
ther magnetic or superconductive 0, U S Si ■ Ear- 
lier, a concept of fermion condensation quantum phase 
transition (FCQPT) preserving quasiparticles and inti- 
mately related to the unlimited growth of M*, had been 
suggested 0, II]. Further studies show that it is capa- 
ble to deliver an adequate theoretical explanation of vast 
majority of experimental results in different HF metals 
EH- In contrast to the Landau paradigm based 



on the assumption that M* is a constant, in FCQPT ap- 
proach M* strongly depends on T, x, B etc. Therefore, 
in accord with numerous experimental facts the extended 
quasiparticles paradigm is to be introduced. The main 
point here is that the well-defined quasiparticles deter- 
mine as before the thermodynamic and transport prop- 
erties of strongly correlated Fermi-systems, while M* be- 
comes a function of T, x, B etc [lllTll, El [H. The 
FCQPT approach had been already successfully applied 
to describe the thermodynamic properties of such differ- 
ent strongly correlated systems as '^He on one side and 
complicated heavy-fermion (HF) compounds on the other 

side [nilillil[illil,[ii. 

In this letter, we analyze the non-Fermi liquid behav- 
ior of strongly correlated Fermi systems us ing the density 
functional theory of fermion condensation . We derive 
equations for the effective mass of quasiparticles in both 
homogeneous systems and solids with a lattice, and show 
that extended quasiparticles paradigm is valid, while the 
dependence of the effective mass on T, x, B etc gives rise 
to the NFL behavior. The obtained results are illustrated 
with calculations of the thermodynamic and transport 
functions of strongly correlated Fermi-systems. Possi- 
ble energy scales in these functions are discussed. We 
demonstrate that our calculations of the heat capacity 
C/T, magnetization M, errergy scales, longitudinal mag- 
netoresistance (LMR) and magnetic entropy S{B) are in 
good agreement with strikiiig recent facts collected on 
the HF metal YbRhaSiz [11, 111 HI. 



II. EQUATION FOR THE EFFECTIVE MASS 
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At first, consider HF liquid at T = characterized 
by the effective mass M*. Upon applying well-known 
Landau equation, we can relate M* to the bare electron 
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Here is the density of states of a free electron gas, M 
is the bare mass, x = p^/Sn"^ is the number density, pp 
is the Fermi momentum, and (x) is the p-wave compo- 
nent of Landau interaction amplitude F. When at some 
critical point x = Xc, F^{x) achieves certain threshold 
value, the denominator in Eq. ^ tends to zero so that 
the effective mass diverges at T = [2l|, and the system 
undergoes FCQPT. It follows from Eq. ^ that beyond 
the critical point Xc, the effective mass becomes nega- 
tive. To avoid an unstable and physically meaningless 
state with a negative effective mass, the system must un- 
dergo a quantum phase transition at the quantum critical 
point X = Xc, which is FCQPT The asymmetrical 

phase behind the quantum critical point is determined by 
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here E is the ground state energy, /i is a chemical po- 
tential, and n{p) is the occupation numbers of quasipar- 
ticles. The main result of such reconstruction is that 
instead of Fermi step, we have < n{p) < 1 in certain 
range of momenta pi < p < Pf. Accordingly, the single- 
particle spectrum 
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in the above momenta interval becomes flat, e(p) = fi, 
and this state is known as a fermi condensate (FC) state 
Due to the above peculiarities of the n(p) function, 
FC state is characterized by the superconducting order 
parameter k(p) = ^n(p)(l — n{p)). 

To derive equation determining the effective mass, we 
employ the density functional theory of superconducting 
state [22. In that case, the ground state energy E be- 
comes the functional of the occupation numbers and the 
function of the number density x, E = E\n( p),x \, while 
Eq. ([31) gives the single-particle spectrum Upon 
differentiating the both sides of Eq. ^ with respect to 
p and after some algebra and integration by parts, we 
obtain 
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Here, F{p, pi) = 5"^ E / 5n{p)5n{pi) is the Landau ampli- 
tude. To calculate the derivative de[p)/dp, we employ 
the functional representation 
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It is seen directly from Eq. Q that the effective mass is 
given by the well-known Landau equation 
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For simplicity, we ignore the spin dependencies. To cal- 
culate M* as a function of T, we construct the free energy 
F — E — TS, where the entropy S is given by 



5 = -2 / Hp) ln(n(p)) + (1 - n(p)) ln(l - n(p))] 



dp 

(2^' 
(7) 

which follows from general combinatorial reasoning []|. 
Minimizing F with respect to n{p), we arrive at the 
Fermi-Dirac distribution, 
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Due to the above derivation, we conclude that Eqs. @ 
and © are exact ones and allow us to calculate the be- 
havior of both d£{p)/dp and M* in the vicinity of FC- 
QPT where the well-defined quasiparticles determine the 
low-tcmperature physics, while M* becomes a divergent 
function of T, B and a; As we wiU see 

this feature of M* forms the NFL behavior observed in 
measurements on HF metals. 



III. SCALING BEHAVIOR OF THE EFFECTIVE 
MASS 

To avoid difficulties associated with the anisotropy 
generated by the crystal lattice of solids, we study the 
universal behavior of heavy-fermion metals using the 
model of the homogeneous heavy-electron (fermion) liq- 
uid. The model is quite meaningful because we consider 
the universal behavior exhibited by these materials at 
low temperatures, a behavior related to power-law diver- 
gences of quantities such as the effective mass, the heat 
capacity, the magnetization, etc. These divergences and 
scaling behavior of the effective mass, or the critical expo- 
nents that characterize them, are determined by energy 
and momentum transfers that are small compared to the 
Debye characteristic temperature and momenta of the 
order of the reciprocal lattice cell length a~^. Therefore 
quasiparticles are influenced by the crystal lattice aver- 
aged over big distances compared to the length a. Thus, 
we can substitute the well-known jelly model for the lat- 
tice as it is usually done, for example, in the fluctuation 
theory of second order phase transitions. 

The schematic phase diagram of HF liquid is reported 
in Fig. [TJ Magnetic field B is taken as the control pa- 
rameter. In fact, the control parameter can be pressure 
P or doping (the number density) x etc as well. At 
B = Bco, FC takes place leading to a strongly degen- 
erated state, where Bco is a critical magnetic field, such 
that at -B > Bco the system is driven towards its Lan- 
dau Fermi liquid (LFL) regime. In our simple model Bcq 
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is a parameter. The FC state is captured by the super- 
conducting (SC), ferromagnetic (FM), antiferromagnetic 
(AFM) etc. states lifting the degeneracy [13, [HI- Be- 
low we consider the HF metal YbRh2Si2. In that case, 
Bco ^ 0.06 T {B±c) and at T = and B < B^o the 
AFM state takes place [20|. At elevated temperatures 
and fixed magnetic field the NFL regime occurs, while 
rising B again drives the system from NFL region to LFL 
one as shown by the dash-dot horizontal arrow in Fig. [1] 
Below we consider the transition region when the system 
moves from NFL regime to LFL one along the horizon- 
tal arrow and it moves from LFL regime to NFL one 
along the vertical arrow as shonw in Fig. [1] The inset to 
Fig. [T] demonstrates the behavior of the normalized effec- 
tive mass M'^ = M* /M^j versus normalized temperature 
T/v = T/Tm, where M^^ is the maximum value that M* 
reaches at T = Tm- The T^"^!"^ regime is marked as NFL 
one since the effective mass depends strongly on tem- 
perature. The temperature region T ~ Tm signifies the 
crossover between the LFL regime with almost constant 
effective mass and NFL behavior, given by T"^/^ depen- 
dence. Thus temperatures T ~ Tm can be regarded as 
the crossover region between LFL and NFL regimes. 




Control parameter, magnetic field B 

FIG. 1: Schematic phase diagram of HF metals. BcO 
is magnetic field at which the effective mass divergences. 
SC, FM, AFM denote the superconducting (SC), ferromag- 
netic (FM) and antiferromagnetic (AFM) states, respectively. 
At B < Bco the system can be in SC, FM or AFM states. The 
vertical arrow shows the transition from the LFL regime to 
the NFL one at fixed B along T with M* depending on T . 
The dash-dot horizontal arrow illustrates the system moving 
from NFL regime to LFL one along B at fixed T . The in- 
set shows a schematic plot of the normalized effective mass 
versus the normalized temperature. Transition regime, where 
Mj^ reaches its maximum value M\i at T = Tm, is shown 
by the hatched area both in the main panel and in the inset. 
The arrows mark the position of inflection point in Mj^ and 
the transition region. 

To explore a scaling behavior of A/*, we write the 
quasiparticle distribution function as ?t.i(p) — n{p,T) — 
7T,(p), with n{p) is the step function, and Eq. ^ then 
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At FCQPT the effective mass M* diverges and Eq. dH) 
becomes homogeneous determining M* as a function of 
temperature 
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while the system exhibits the NFL behavior [10|, |l3[ . If 
the system is located before FCQPT, M* is finite, at 
low temperatures the system demonstrates the LFL be- 
havior that is M*(T) ~ M* + aiT^, with oi is a con- 
stant, see the inset to Fig. [T] Obviously, the LFL behav- 
ior takes place when the second term on the right hand 
side of Eq. ^ is small in comparison with the first one. 
Then, at rising temperatures the system enters the tran- 
sition regime: M* grows, reaching its maximum M^j at 
T — Tm, with subsequent diminishing. Near tempera- 
tures T > Tm the last "traces" of LFL regime disappear, 
the second term starts to dominate, and again Eq. ([9]) 
becomes homogeneous, and the NFL behavior restores, 
manifesting itself in decreasing of M* as T~'^/^. When 
the system is near FCQPT, it turns out that the solution 
of Eq. dni) M*{T) can be well approximated by a simple 
universal interpolating function [lol . [isl . fisj . The interpo- 
lation occurs between the LFL (M* ~ M* aiT^) and 
NFL (M* oc T-2/3) regimes thus describing the above 
crossover [l3, [l^- Introducing the dimensionless vari- 
able y = Tat = T /Tm, we obtain the desired expression 
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Here M'^ = M* /Ml^ is the normalized effective mass, 
Co = (1 + C2)/(l + ci), ci and C2 are fitting parameters, 
parameterizing the Landau amplitude. 

It is possible to transport Eg. ^ to the case of the 
application of magnetic fields |13l . [l5| . The application 
of magnetic field restores the LFL behavior so that 
depends on B as 
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Tm oc ^jlb{B - Bco), 
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where /is is the Bohr magneton [10|, |l3|, [l5|. Employ- 
ing Eqs. HI]) and ((T3]) to calculate Mlj and Tm, we 
conclude that Eq. ((TT|) is valid to describe the normal- 
ized effective mass in external fixed magnetic fields with 
y = T/{B - Bco)- On the other hand, Eq. ^ is valid 
when the applied magnetic field becomes a variable, while 
temperature is fixed T ~ Tf. In that case, as seen from 
Eqs. (fTO|) . (fTT|) and p^ . it is convenient to rewrite both 
the variable as y — {B — Bco)/Tf, and Eq. p3|) as 



fiB^BM - Bco) oc Tf. 
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It follows from Eq. pT|) that in contrast to the Landau 
paradigm of quasiparticles the effective mass strongly de- 
pends on T and B. As we will see it is this dependence 
that forms the NFL behavior. It follows also from Eq. 
(fTTj) that a scaling behavior of M* near FCQPT point is 
determined by the absence of appropriate external physi- 
cal scales to measure the effective mass and temperature. 
At fixed magnetic fields, the characteristic scales of tem- 
perature and of the function M* (T, B) are defined by 
both Tm and Mlj respectively. At fixed temperatures, 
the characteristic scales are {Bm — Bco) and M^j. It 
follows from Eqs. and that at fixed magnetic 

fields, Tm 0, and M^j oo, and the width of the 
transition region shrinks to zero as B Bcq when these 
are measured in the external scales. In the same way, it 
follows from Eqs. pil)) and (|14p that at fixed tempera- 

cxD, and the width of 



0, and M 
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tures, {Bm - B^ 
the transition region shrinks to zero as T/ ^0. 

A few remarks are in order here. As we shall see in 
Subsections IV C and D, in some cases temperature or 
magnetic field dependencies of the effective mass or of 
other observable like the longitudinal magnetoresistance 
do not have "peculiar points" like maximum. The nor- 
malization are to be performed in the other points like 
the inflection point shown in the inset to Fig. [TJ Such 
a normalization is possible since it is established on the 
internal scales. 



IV. NON-FERMI LIQUID BEHAVIOR IN 

YbRh2Si2 

In what follows, we compute the effective mass and 
employ Eq. (fTTj) for estimations of considered values. To 
obtain the effective mass M* (T, B), we solve Eqs. (HI and 
(|9]) with special form of Landau interaction amplitude, 
see Refs. 0, for details. Choice of the amplitude is 
dictated by the fact that the system has to be in the FC- 
QPT point, which means that first two p-derivatives of 
the single-particle spectrum e(p) should equal zero. Since 
first derivative is proportional to the reciprocal quasipar- 
ticle effective mass I /A/*, its zero (where 1/M* = and 
the effective mass diverges) just signifies FCQPT. Zeros 
of two subsequent derivatives mean that the spectrum 
e(p) has an inflection point at Fermi momentum pp so 
that the lowest term of its Taylor expansion is propor- 
tional to {p - pf)^ [il- After solution of Eq. the 
obtained spectrum had been used to calculate the en- 
tropy S{B,T), which, in turn, had been recalculated to 
the effective mass M* (T, B) by virtue of well-known LFL 
relation M*{B,T) = S{B,T)/T. We note that our cal- 
culations conflrm the validity of Eq. (fTT|) . 



A. Heat capacity and the Sommerfeld coefficient 

Exciting measurements of C/T oc M* on samples of 
the new generation of YbRh2Si2 in different magnetic 



fields _B up to 1.5 T [l8| allow us to identify the scaling 
behavior of the effective mass M* and observe the dif- 
ferent regimes of M* behavior such as the LFL regime, 
transition region from LFL to NFL regimes, and the NFL 
regime itself. A maximum structure in C/T oc A/|j at 
Tm appears under the application of magnetic field B 
and Tm shifts to higher T as B is increased. The value 
of C jT = 7o is saturated towards lower temperatures 
decreasing at elevated magnetic field, where 70 is the 
Sommerfeld coefficient [l8| . 

It follows from Section III, that the transition region 
corresponds to the temperatures where the vertical arrow 
in the main panel of Fig. [1] crosses the hatched area. 
The width of the region, being proportional to Tm oc 
(B — Bcq) shrinks, Tm moves to zero temperature and 
7o oc M* increases as i? — > Bco- These observations are 
in accord with the facts UM- 
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FIG. 2: The normalized effective mass M'^ extracted from 
the measurements of the specific heat C/T on YbRh2Si2 in 
magnetic fields B shown in the left down corner [l^. Our 
calculations are depicted by the solid curve tracing the scaling 
behavior of M^. 



To obtain the normalized effective mass M , the max- 
imum structure in C/T was used to normalize C/T, and 
T was normalized by Tm- In Fig. [5]the obtained M'^ as a 
function of normalized temperature T/v is shown by geo- 
metrical figures, our calculations carried out as described 
above are shown by the solid line. Figure [2] reveals the 
scaling behavior of the normalized experimental curves 
- the curves at different magnetic fields B merge into a 
single one in terms of the normalized variable y = T /Tm- 
As seen from Fig. [2l the normalized mass extracted 
from the measurements is not a constant, as would be for 
a LFL, and shows the scaling behavior given by Eq. pdjl 
over three decades in normalized temperature. The two 
regimes (the LFL regime and NFL one) separated by the 
transition region, as depicted by the hatched area in the 
inset to Fig. [1] are clearly seen in Fig. [2] illuminating 
good agreement between the theory and facts. 
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B. "Average" magnetization 




1 2 3 4 5 



Normalized magnetic field 



It is seen from Fig. [3] that our calculations are in good 
agreement with the facts, and all the data exhibit the 
kink (shown by arrow) at Bn ~ 1 taking place as soon 
as the system enters the transition region corresponding 
to the magnetic fields where the horizontal dash-dot ar- 
row in the main panel of Fig. [T] crosses the hatched area. 
Indeed, as seen from Fig. [31 at lower magnetic fields M 
is a linear function of B since M* is approximately inde- 
pendent oi B. It follows from Eq. (|12p that at elevated 
magnetic fields M* becomes a diminishing function of B 
and generates the kink in M{B) separating the energy 
scales discovered in Ref. [l^. Then, it seen from Eq. 

that the magnetic field Bk — Bm at which the kink 
appears shifts to lower B as Tf is decreased. 

C. Longitudinal magnetoresistance 



FIG. 3: The field dependence of the normalized "average" 
magnetization M = M + Bx is shown by squares and ex- 
tracted from measurements collected on YbRu2Si2 [3. The 
kink (shown by the arrow) is clearly seen at the normalized 
field Bn = B/Bk — 1. The solid curve and stars (see text) 
represent our calculations. 

Consider now an "average" magnetization M = Bx + 
M as a function of magnetic field B at fixed temperature 
T = Tf, where x is the magnetic susceptibility and M is 
the magnetization. 



(15) 



M{B,T)^ / x{b.T)dh, 
Jo 

where the magnetic susceptibility is given by [l| 

(3M*{B,T) 



xiB,T) = 



1 + 



(16) 



Here, /3 is a constant and Fq is the Landau amplitude 
related to the exchange interaction pj. In the case of 



strongly correlated systems Fq > —0.9 [2l|. Therefore, 
as seen from Eq. (jl6[) . due to the normalization the co- 
efficients P and (1 + Fq) drops out from the result, and 
X oc M*. To obtain M, we calculate M* as a function 
of B at fixed Tf . The obtained curves of M exhibit en- 
ergy scales separated by kinks aX B = Bk- As seen from 
Fig. O the kink is a crossover point from the fast to slow 
growth of M at rising magnetic field. We use Bk and 
M{Bk) to normalize B and M respectively. 

The normalized M vs the normalized field -Bjv = 
B /Bk are shown in Fig. [31 Our calculations are depicted 
by the solid line. The stars trace out our calculations of 
M with M*{y) extracted from the data C/T shown in 
Fig. [51 The calculation procedure deserves a remark 
here. Namely, in that case M* depends on y = T/Tm 
with Tm is given by Eq. (|13p. On the other hand, we 
can consider y — (B — Bco)/Tf as it is shown in Section 
III, and take the data C/T as a function of y. 
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FIG. 4: Magnetic field dependence of the normalized (in the 
infiection point shown by the arrow, see text for details) mag- 
netoresistance R'^ versus normalized magnetic field, ii^ was 
extracted from LMR of YbRh2Si2 at difi'erent temperatures 
[l9l | listed in the legend. The solid line represents our calcu- 
lations. 

Consider a longitudinal magnetoresistance (LMR) 
p{B,T) = po + AT^ as a function of B at fixed Tf. 
In that case, the classical contribution to LMR due to 
orbital motion of carriers induced by the Lorentz force 
is small, while the Kadowaki- Woods relation [H, [23 |. 
K = A/^Q (X A/x^ = const, allows us to employ M* 
to construct the coefficient A [2^, since 70 oc % a M* . 
As a resuh, p{B,T) - po oc {M*f. Fig. H reports the 
normalized magnetoresistance 



pjy) - Po 

Pinf 
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vs normalized magnetic field y — B/Binf at different 
temperatures, shown in the legend. Here pinf and Bm/ 
are LMR and magnetic field respectively taken at the 
inflection point marked by the arrow in Fig. [4l 

The normalization procedure deserves a remark here. 
Namely, since the magnetic field dependence of both the 
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calculated M* and LMR does not have "peculiar points" 
like maximums, the normalization have been performed 
in the corresponding inflection points. To determine the 
inflection point precisely, we first differentiate p(B,T) 
over B, find the oxtremum of derivative and normalize 
the values of the function and the argument by their val- 
ues in the inflection point. Then, both theoretical (shown 
by the solid line) and experimental (marked by the ge- 
ometrical figures) curves have been normalized by their 
inflection points, which also reveal the universal behavior 
- the curves at different temperatures merge into single 
one in terms of the scaled variable y and show the scaling 
behavior over three decades in the normalized magnetic 
field. 

The transition region at which LMR starts to decrease 
is shown in the inset to Fig. [T]by the hatched area. Obvi- 
ously, as seen from Eq. (|14p , the width of the transition 
region being proportional to Bm — Binf decreases as the 
temperature Tf is lowered. In the same way, the inflec- 
tion point of LMR, generated by the inflection point of 
M* shown in the inset to Fig. [T]by the arrow, shifts to 
lower B as Tf is decreased. All these observations are in 
excellent agreement with the facts [l9l |. 



D. Magnetic entropy 

The evolution of the derivative of magnetic entropy 
dS{B,T)/dB as a function of magnetic field B at fixed 
temperature Tf is of great importance since it allows us 
to study the scaling behavior of the derivative of the ef- 
fective mass TdM*{B,T)/dB oc dS{B,T)/dB. While 
the scaling properties of the effective mass M*{B, T) can 
be analyzed via MLR as we have shown in Subsection C. 
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NORMALIZED MAGNETIC FIELD 



FIG. 5: Normalized magnetization difference divided by tem- 
perature increment {AM / AT)n vs normalized magnetic field 
at fixed temperatures (listed in the legend in the upper left 
corner) is extracted from the facts collected on YbRh2Si2 [20i |. 
Our calculations of the normalized derivative {dS/dB)N — 
(AM/Ar)jv vs normalized magnetic field are given at fixed 
dimensionless temperatures T/n (listed in the legend in the 
upper right corner). All the data are shown in the geometrical 
figures depicted in the legends. 



As seen from the consideration in Subsection C and 
from Eqs. pT|) and at y < 1 the derivative 

-dM]si{y)/dy oc y with y = {B - Bco)/{Binf - Bco) oc 
{B — Bco)/Tf. We recall that the effective mass as a 
function of B does not have the maximum, see Sub- 
section C. At elevated y the derivative —dMM{y)/dy 
possesses a maximum at the inflection point and then 
becomes a diminishing function of y. Upon using the 
variable y = {B ~ BcQ)/Tf, we conclude that at de- 
creasing temperatures, the leading edge of the function 
~dS/dB oc ~TdM* / dB becomes steeper and its maxi- 
mum at {Binf — Bcq) oc Tf is higher. These observations 
are in quantitative agreement with striking measure- 
ments of the magnetization difference divided by tem- 
perature increment, — AM/ AT, as a function of magnetic 
field at fixed temperatures Tf collected on YbRh2Si2 [20j . 
We note that according to the well-know thermodynamic 
equality dM/dT = dS/dB, and AM/ AT ~ dS/dB. 

To carry out a quantitative analysis of the scaling 
behavior of —dM*{B,T)/dB, we calculate as described 
above the entropy S{B,T) as a function of B at fixed 
dimensionless temperatures T/Z/i shown in the upper 
right corner of Fig. [51 This figure reports the normal- 
ized (dS/dB) AT as a function of the normalized magnetic 
field. The function (dS/dB) jv is obtained by normal- 
izing (—dS/dB) by its maximum taking place at Bm, 
and the field B is scaled by Bm- The measurements 
of —AM/ AT are normalized in the same way and de- 
picted in Fig. O as {AM/ AT)jv versus normalized field. 
It is seen from Fig. [5] that our calculations are in excel- 
lent agreement with the facts and both the experimen- 
tal functions {AM/AT)n and the calculated {dS/dB)^ 
ones show the scaling behavior over three decades in the 
normalized magnetic field. 



V. SUMMARY 

We have analyzed the non-Fermi liquid behavior of 
strongly correlated Fermi systems using the density func- 
tional theory of fermion condensation and derived equa- 
tions for the effective mass of quasiparticles in both ho- 
mogeneous systems and solids with a lattice, and showed 
that extended quasiparticles paradigm is strongly valid, 
while the dependence of the effective mass on temper- 
ature, number density, applied magnetic fields etc gives 
rise to the NFL behavior. The obtained results are illus- 
trated with calculations of the thermodynamic and trans- 
port functions of strongly correlated Fermi-systems. Pos- 
sible energy scales in these functions are discussed. We 
have demonstrated that our comprehensive theoretical 
study of the heat capacity, magnetization, energy scales, 
the longitudinal magnetoresistance and magnetic entropy 
are in good agreement with the outstanding recent facts 
collected on the HF metal YbRh2Si2. 
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